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Abstract. We employ a new method for studying compositionally disordered ferroelectric ox- 
ides. This method is based on the virtual crystal approximation (VCA), in which two or more 
component potentials are averaged into a composite atomic potential. In our method, we con- 
struct a virtual atom with the correctly averaged atomic size and atomic eigenvalues. We have 
used our new method to study the composition dependent phase transition in Pb(Zri_ x Ti 2 .)03 
lying between x = 0.5 and x — 0.4. We correctly predict the experimentally determined 
phase transition from the tetragonal phase to a low-temperature rhombohedral phase between 
these two compositions. 



INTRODUCTION 



Our recent work [1] has demonstrated the utility of the virtual crystal approxima- 
tion (VCA) [2,3] to study compositionally disordered materials. Following our initial 
VCA study of a stress-induced phase transition in Pb(Zr 1 _ a .Ti a; )03, we and other authors 
have extended the VCA to studies of compositional phase transitions [4], temperature- 
dependent phase transitions [5], and dynamical properties. [6] Prior to this the VCA has 
had some success in providing qualitative, and in some instances quantitative agreement 
with large-scale solid-state calculations. 

Previous implementations of the VCA have focused on the averaging of two or more 
component potentials in separable form at the solid-state level. [7,8] By averaging the 
potentials, it is difficult to assess the quality of the resulting VCA potential. In some 
studies, these traditional VCA potentials have yielded unphysical results. [9,10] Spec- 
ulation has been made in the literature that it is the inability of these traditional VCA 
potentials to capture the important differences in chemical bonding and ionicity accu- 
rately that has given the poor results. [9] 



METHODOLOGY 



Due to the inconsistencies in the literature regarding the quality of VCA results, we 
have recently formulated a new approach for the construction of VCA potentials. Our 
method represents a simple and intuitive way to construct VCA potentials and assess 
their quality. It is a severe departure from the more traditional methods for VCA potential 
construction. We have previously presented a thorough description of our method for the 
construction of VCA atoms . [ 1 ,4] For brevity, we will only present here the salient points 
of our method. 

• Our method averages the component atoms at the all-electron level, such that all- 
electron eigenvalues, potentials, wave functions and charge density are all com- 
puted for the virtual atom. Since the solution of the all-electron atom is common 
to all pseudopotential procedures, our method is general enough to be used in all 
types of pseudopotential construction. 

Operationally we find the all-electron wave functions and eigenvalues for the atoms 
we want to average. We then determine the properly averaged eigenvalues for the 
VCA atom valence states from the eigenvalues of the all-electron atoms. We then 
construct a bare Coulombic potential from the properly averaged bare nuclear po- 
tentials of the component atoms and a properly averaged core charge density. In 
doing so, we have insured that the resulting VCA potential will have the correct av- 
eraged size. Using this averaged nuclear potential and frozen-core charge density, 
we construct new self-consistent wave functions for the VCA atom valence states. 
We require that these new wave functions solve the Kohn-Sham equations in the 
valence region and give the properly averaged eigenvalues. 

• From this new set of potentials and wave functions, we construct an optimized 
pseudopotential. [11] This method enforces the exact agreement of the VCA eigen- 
values with the averaged all-electron eigenvalues for the reference state only. For 
any other electronic configuration, we have made no requirements on the quality 
of the potential. Recently, we have formulated the designed nonlocal pseudopo- 
tential method, [12] in which we exploit the inherent arbitrary separation of the 
local potential and semilocal correction terms required for the Kleinman-Bylander 
separable form. [13] By adjusting the form of the local potential and therefore the 
correction terms, we may dramatically improve the transferability of the potential 
without affecting the exact agreement at the reference configuration. This transfer- 
ability improvement over a wide range of electronic configurations insures correct 
ionicity and polarizability of the virtual atom. 

It is our assertion that the preservation of the proper averaged size of the VCA po- 
tential and a high level of transferability will provide accurate VCA atoms which yield 
excellent agreement with superlattice calculations. We have previously demonstrated 
the effectiveness of our new VCA approach. [1,4] The remainder of this paper will fo- 
cus on the construction and use of two VCA potentials to study the compositional phase 
transition in Pb(Zr 1 _ 2 ,Ti 2 ,)0 3 . 



RESULTS AND DISCUSSION 



We have applied our new VCA method to the 7ii- x Ti x virtual atom for x = 0.5 and 
x — 0.4. All atomic energy calculations were done with the local density approximation 
and optimized [11] and designed nonlocal pseudopotential [12] methods were employed. 
The generation parameters for the two VCA potentials are included in Table 1 . For all 
atoms, semi-core states were included as valence. It is important to note that although we 
have included multiple s-channel states, only one s nonlocal projector is used. For both 
VCA atoms, we have used the s-potential with the addition of one or two square-step 
potentials as the local potential. 

The transferability testing results for the two VCA potentials are contained in Table 2. 
From the table, it is evident that our new VCA pseudopotentials are highly transferable 
at either value of x. 

We use both VCA potentials in solid-state calculations for three structural phases of 
Pb(Zr 1 _ 2 .Ti 2 .)0 3 . The electronic wave functions are expanded in a plane- wave basis us- 
ing a cutoff energy of 50 Ry. We included the 5d shell as valence for the Pb atom as 
well as including scalar relativistic effects. Brillouin zone integrations are approximated 
accurately as sums on a 4x4x4 Monkhorst-Pack fc-point mesh. [14] 

We have included three structurally distinct phases of Pb(Zri_ x Ti :r )03 in the current 
studies of compositions between x = 0.5 and x = 0.4: a tetragonal Ti-rich Pimm and 
two rhombohedral Zr-rich (R3c and R3m) phases. [15,16] Around 400-500K, the rhom- 
bohedral region exhibits a boundary between R3c (low-temperature) and R3m (high- 
temperature) phases, which depends weakly on composition. [17,18] The R3c phase 
shows complex oxygen octahedral tilting, which doubles the primitive unit cell to ten 
atoms. [19] 

We have completed full electronic and structural relaxations for five-atom unit cells for 
the tetragonal and high-temperature rhombohedral phases. Due to the complex oxygen 
octahedral distortions, we have used a 10-atom unit cell for the low-temperature rhom- 



TABLE 1. Construction parameters for the Zri-^Tia; virtual crys- 
tal (VCA) designed nonlocal pseudopotentials. The VCA poten- 
tials were generated with the methods described in references [11], 
[12] and the method described in text. Core radii (r c ) are in atomic 
units, q c parameters are in Ry 1 / 2 , step widths are in atomic units 
and step heights are in Ry. 



Atom 



Step 
Range 



Step 
Height 



VCA (s 2 p 6 d°) 1.38,1.51,1.40 

x = 0.5 



7.07 0.00-0.56 
0.56-0.79 



1.18 
1.34 



VCA(sVd°) 1.27,1.38,1.61 7.07 0.00-1.25 10.32 

x = 0.4 



TABLE 2. Configuration testing for Zri-^Tia; 
virtual crystal (VCA) atoms generated with the 
method described in text. Averaged eigenvalues 
and total energy differences are given for Zr and Ti 
pseudopotentials (PS) for two values of x. The PS 
construction parameters have been presented else- 
where. [1,4] Errors are computed as a difference 
between VCA and averaged PS results. All ener- 
gies are in Ry. 
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bohedral phase. For all rhombohedral calculations, we have neglected the small shear 
relaxations. 

In Figure 1, we show the equations of state for the three phases of Pb(Zri_a;Ti x )03 at 



x = 0.5. We find the tetragonal phase is the ground-state structure for this composition 
lying approximately 0.05 eV lower in energy than both rhombohedral phases. This is 
correct energy ordering found experimentally. [15,16] 
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FIGURE 1. Equations of state for the tetragonal (solid line), high-temperature rhombohedral (dot- 
ted-dashed line) phase and low-temperature rhombohedral (dotted line) phases of Pb(Zri_ x Ti 2; )03 for 
x = 0.5. The heights and energies are for a 40-atom unit cell. 

In Figure 2, we show the equations of state for the three phases of PbCZr^-j. Ti rr )0 3 
at x = 0.4. We find both rhombohedral phases have shifted down in energy relative to 
the tetragonal phase. In fact, the low-temperature rhombohedral phase has shifted by 
approximately 0.09 eV in energy and therefore becomes the ground-state structure for 
this composition of PbCZri.^Ti^OOs. This result is in direct agreement with the experi- 
mental findings [15,16] and demonstrates the ability of our new VCA method to locate 
and predict compositional phase transitions. 

CONCLUSIONS 

In this paper, we have applied our new method for constructing a virtual crystal pseu- 
dopotential to the Zr and Ti atoms. Potentials constructed with our method not only 
possess the properly averaged atomic size but also a high level of accuracy in describing 
energetic differences due to changes in both ionicity and polarizability. This method is 
based on averaging all-electron information and therefore makes the method applicable 
to all types of pseudopotential construction algorithms. We have used our new method 
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FIGURE 2. Equations of state for the tetragonal (solid line), high-temperature rhombohedral (dot- 
ted-dashed line) phase and low-temperature rhombohedral (dotted line) phases of Pb(Zri_ K Ti a ;)03 for 
x = 0.4. The heights and energies are for a 40-atom unit cell. 

to construct two Zri_ x Ti x VCA potentials at x = 0.5 and x = 0.4. With these potentials, 
we have studied three structural phase of Pb(Zri ^ Ti x )03, finding that the tetragonal 
phase is the ground-state structure for Pb(Zri_ :E Ti :E )03 at x = 0.5. When moving to 
x = 0.4, we find that the low-temperature rhombohedral phase becomes the ground- 
state. This finding is in excellent agreement with experiment and represents the first ab 
initio determination of a compositionally-dependent phase transition in a ferroelectric 
oxide. 

Recently yet another phase (monoclinic Cm) has been experimentally found at low 
temperatures for compositions between x = 0.45 and x = 0.5. [20] We plan to extend 
the current research by studying the monoclinic phases with our VCA atoms. 
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